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Abstract 

From the analysis of Supernova la data alongwith Observational Hubble Data (OHD) and 
Baryon Acoustic Oscillation (BAO) data, we attempt to find out the nature of a scalar potential 
that may be responsible for the Dark Energy of the universe. We demonstrate that in order to 
explain the varying dark energy equation of state (u>x(z)) as obtained in a model independent 
way from the analyses of observational data , we need to invoke a quintom scalar field having 
a "quintessence" part for ujx(z) > — 1 and a "phantom" part for ujx(z) < — 1. We consider a 
Gaussian type potential for these scalar fields and compare the dark energy equation of state 
derived from such potential with the one computed from the data analysis. 
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1 Introduction 



The recent acceleration of the expansion rate of the universe that has been inferred from the 
observation and luminosity distance measurements of different Supernova la (SNe la) [1] at different 
redshifts (z) is thought to be caused by the existence of a mysterious energy called dark energy. 
Different other observational evidence (such as WMAP [2], Cosmic Microwave Background (CMB) 
[3]) etc. suggest that this dark energy accounts for almost 73% of the content of the universe. 
Dark energy is a hypothetical energy component with a negative pressure responsible for causing 
the universe to experience this accelerated expansion. Although there are efforts to understand its 
nature and origin, dark energy still remains an enigma. 

An attempt to explain dark energy involvs considering a slowly varying potential V{4>) of a 
scalar field (f>. Hence for such considerations the knowledge of <f> and V(4>) are necessary in order 
to get an insight of dark energy and its equation of state. In FRW cosmology, it can be shown 
that in cosmological constant scenario , the dark energy equation of state remains constant at the 
value -1 all along the evolution history of the universe. However if dark energy varies with time 
or equivalently with redshift z, a possible explanation may be given by a scalar field (quintessence 
field) for dark energy equation of state u)(z) > —1. But one needs to invoke a phantom scalar field 
a in case w(z) falls below —1. 

In the present work, we have considered five different SNe la data sets alongwith Observational 
Hubble Data (OHD) [4] and Baryon Acoustic Oscillation Data (BAO) [5]. The data sets are made 
available from Riess et al (2007) [6], Wood-Vasey et al (2007) [7], Davis et al[8], Kowalski et al 
(2008) [9], Kessler et al (2009) [10] and Amanullah et al (2010) (also known as Union2 data) [11]. 
The maximum span of redshift z for these compilations for distance modulus fi(z) (related to the 
luminosity distance di(z)) are in the range 0.001 < z < 1.76. A parametric form for <1l{z) is chosen 
whose parameters are obtained from a x 2 minimization of various observational data mentioned 
above. The dark energy equation of state can be derived analytically using an analytic form of 
di(z), such as the present parametric form. With the fitted values of the parameters the dark energy 
equation of state cux(z) can therefore be calculated. On the otherhand, since cIl(z) is related to 
the Hubble parameter H (z) by an analytical expression one can use the same parametric form for 
cIl(z) (with parameters from \ 2 analyses) and evaluate both the quintessence scalar field (j> and the 
phantom field a and their variations with the redshift z (see later) . With chosen form of potentials 
for both (j> and a, one can then compute the dark energy equation of state (uj(z)) from the 4> an d 
a values obtained from the analysis. In this work, we have taken the form of this potential as 
V = V{(f>) + V{a), where V(<p) and V(a) are chosen to be V(<f>) = V exp[Ai(l + (<f>(z) - (p ) 2 /M^)} 
and V(a) = Voexp[A2(l + (cr(,z) — <to) 2 /M5)], where 0o(o"o) is the value of 4>{cr) at present epoch and 
M p i is the Planck mass. We refer the dark energy equation of state computed from the scalar fields 
with this potential, as oj(z). The uj(z) thus obtained is then compared with the ojx(z) computed 
directly from the analysis of ^l{z) data with the assumption that the present universe is spatially 
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flat and contains only matter and dark energy. 

The paper is organised as follows. In Sect. 2 we discuss the dark energy equation of state as 
derived from FRW cosmology. In Sect. 3 we propose a quintom scalar formalism of dark energy. 
The calculational procedures are given in Sect. 4. Finally in Sect. 5 we furnish our calculational 
results and discussions. 

2 Brief overview of model independent reconstruction of dark en- 
ergy equation of state 

For a homogeneous and isotropic universe the Friedmann-Robertson- Walker metric is given by 



ds 2 = -dt 2 + a 2 (t) 



1 - Kr 2 



+ r 2 (d6 2 + sin 2 6d4> 2 ) 



(1) 



where a(t) is the scale factor and K is the curvature parameter. Considering the universe to be a 
perfect fluid characterized by energy density p and pressure p, the Einstein's equation alongwith 
Eq. (1) leads to two independent equations (for spatially flat universe), 

8ttG 



H 2 (t) 
H(t) 



-Pit) 



-4nG(p + p) 

where H(t) = («(*)) is the Hubble parameter that denotes the expansion rate of universe 
From Eq. (2) and Eq. (3) one gets 

alt) 4vrG. 
a(t) 3 

Assuming the universe contains only matter and dark energy, Eq. (2) leads to the result 



(2) 
(3) 



(4) 



H 2 (z) 
H 2 



+ zf + n° x exp ( [ Z 3(1 + u x (z))-^- 
\Jo 1 + z 



(5) 



where the symbol £1 represents energy density or matter density normalised to the critical density 
p c of the universe and thus (&x) represents the matter (dark energy) density parameter at 
present epoch. As mentioned earlier in this section that universe contains only matter and dark 
energy in the redshift range considered here, we have tl m (z) + £lx(z) = 1. From Eq. (5) the 
equation of state of dark energy oux(z) can be derived as, 



-1 + 



2 (l+z) dH(z) _ Q o (1+z) 3 tt2 

1 - o° il+£>i H 2 



(6) 



The SNe la observational data are tabulated as a distance modulus p(z). This is related to 
luminosity distance (1l(z) by 



p(z) 



5log w (d L /Mpc) + 25 



(7) 
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for a redshift value z at which the observation was made. The observed distance modulus n bs( z ) 
is related to the apparent magnitude m ^ s (z) of a SNe la and absolute magnitude M through 
the equation /i bs(-z) = mobs(z) — M. Also for spatially fiat universe Hubble parameter and the 
luminosity distance are related through 



As stated in Sect. 1 a parametric form of cLl{z) is considered and the parameters are obtained 
from a x 2 fit of the observational data. The dark energy equation of state oox(z) can then be 
reconstructed using the Eqns (1 - 5) and Eq. (6). We also mention here that a combined \ 2 fit is 
made with SNe la data, OHD and BAO data. 

3 Dark energy equation of state in a scalar field theory framework 

The acceleration equation (Eq. (4)) suggests that the universe containing a perfect fluid can undergo 
accelerated expansion only if the pressure of the fluid is negative. A suitably chosen real scalar field 
can produce negative pressure in the FRW spacetime and may also lead to time varying equation 
of state of the perfect fluid depending on the potential of the field. Therefore such a real scalar 
field can be a possible candidate for time varying dark energy. The scalar fields naturally arise in 
particle physics theories including string theory. There exists in literature a wide variety of scalar 
field dark energy models. In this work we consider a quintom scalar field model where the two 
scalar fields namely <p(z) and <j(z) are responsible for dark energy and the variation of its equation 
of state with redshift z. Each of the fields has specified form of potential and the two kinds of fields 
are not coupled to each other. We like to show that this simple model can explain the varying dark 
energy equation of state ujx(z) obtained by using the combined \ 2 analyses results of SNe data, 
OHD data, BAO data in Eq. (6). 

3.1 Dark energy and quintessence scalar field 

The action S for Quintessence scalar field is given by (with V (</>), the potential of the field (ft), 



where S m is the matter action. The equation of motion of the spatially homogeneous quintessence 
field is given by 




(8) 



S 




(9) 



4>{t) + 3H(t)<j>(t) + 



dV{(j)) 
dct) 



= 



(10) 



The Eq. (10) will lead to 



1 d 

&irGH{t) dt 



{Hi) 



(11) 
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where #J = S*p (±0 2 (i) + V(<£)) and H 2 = (\<j) 2 {t) + V{<j>) + p m )- The Energy momentum 
tensor for quintessence scalar field can be written as 

T£ = g^dafat-SjiQg^daWrt + Vtt)) (12) 

From Eq. (12) we have, 

P4> = -T ° = ^ 2 + ^) (13) 

P<p = Ti= l -4?-V{<t>) (14) 

Mt) = * = (15) 

where p<p and are respectively the pressure and density due to the scalar field (j). It is evident 
from Eq. (15) that oj^ can vary between -1 and +1. 

3.2 Dark energy and phantom scalar field 



The action for Phantom scalar field a(t) is given by (with V(a), the potential of the field a) 



S m + j (fxyf^g 



l -g^{d^)d v a) + V{a) 



(16) 



where S m is the matter action. The equation of motion of the spatially homogeneous phantom field 
can be written as 

<t(t) + SH(t)&(t) - = 0. (17) 

da 

The above equation leads to 

^ = stoics • < 18 > 

where H 2 = (^cr 2 (t) + V(<r)) and H 2 = ^f- (\a 2 {t) + V(a) + p m ) . The Energy momentum 
tensor for phantom scalar field is given by 

T£ = -grd a <7d v a-Si{^^d a GdpG + V[(T)^ . (19) 

From Eq. (19) we have, 

Pa = -T ° = -V + y(a) (20) 

Pa = Tt = -\a 2 -V(a) (21) 

Pa -^cr^ + K(cr) 
The above equation shows that u a can vary between -1 and -oo. 
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3.3 Quintom scalar field motivated dark energy equation of state 

We consider here a model [12] which contains a negative kinetic scalar field a and a normal scalar 
field 4> with a general potential, described by 



S = S m + J d 4 Xy^ 



(23) 



.16ttG 2" " " 2' 

where S m = J d 4 xC m and C m represents the Lagrangian density of matter fields. Assuming the 
fields to be homogeneous, in a spatially flat FRW cosmological model the equation of motion of the 
fields and the matter density are given by 

i(t) + 3H(t)<j>(t) + ^- = (24) 

a(t) + 3H(t)a{t) - = (25) 

da 

p~ i + 3H{t){p~ i +p 1 ) = 0, (26) 

where the p 7 is the density of fluid (matter) with a barotropic equation of state p 7 = (7 — l)p 7 , 7 
being a constant with < 7 < 2. For radiation, 7 = 4/3 and for dust 7 = 1. The corresponding 
effective pressure and energy densities will be 

V = -\d 2 + l -^-V(a^) (27) 

P = -\ct 2 + \^ + V(<j^) (28) 

In this formalism the Friedmann equation becomes 

H 2 = ^(-\o 2 + V{e) + \^ + V^) + Pl ) 

= Hl + Hl + ^p, (29) 

and the equation of state parameter lo(z) as obtained from the present quintom scalar field formal- 
ism will then be written as 

u, = -O' + f-^fl. (30 ) 
-a 2 + ^ 2 + 2V>,0) 

If the two scalar fields (a and <p) are not directly coupled to each other then V(a, <p) can be written 
as V(a, (f>) = V(a) + V{4>). With this assumption one can readily see that for <j>> a, Eq. (30) gives 
tow> —1 implying the scenario for quintessence scalar field and for <\> < a, oj < — 1 signifying the 
phantom scalar field model. 



4 Calculational Procedure 

The purpose of this work is to show the viability of the present formalism of quintom scalar field 
model in explaining the dark energy and the variation of its equation of state with z. In order 
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to obtain the dark energy equation of state from observational data, a parametric form for di,(z) 
[13, 14] is first considered whose parameters are fixed by the x 2 minimization of the combined data 
(SNe la, OHD, BAO). The matter density at the present epoch Q° m is also made a parameter in 
this x 2 analyses and its value is obtained from the same x 2 minimization. The form of di(z) can 
now be used to calculate Hubble parameter H(z) using Eq. (8). Eq. (6) now readily gives the dark 
energy equation of state wx(z). 

The parametric form of di(z) should respect the conditions that di(z) = at z = and 
di(z) oc z for large z. SNe la data are available for z & 1.76 from various supernova observations. 
A parametrised form of <1l(z) can be written as [13, 14] 



d L {z) = -£r 



z(l + az) 



l + bz 



(31) 



The parameters a, b and Q,% t are obtained by minimising a suitably defined x 2 - In the present case 
the x 2 is defined as [14], 

Xtot = XSN + XBAO + XOHD (32) 

where x!n corresponds to the SNe la data, Xbao defines the x 2 f° r Baryon acoustic oscillation data 
and Xohd represents the x 2 for Observational Hubble Data (OHD). 
x|n [15] is defined as 

X SN (a,&,M ) = 2^ 



In the above, the theoretical distance modulus /x t h(a, b, z) can be obtain from Eq. (7). The above 
equation can be written in the form [15] 

X 2 SN (a,b,M' = B/C) = A-^ r . (33) 



where 



A = ^ ( 51 °gio(^( Q ' 6 '^)) ~ m obs{zi)) 2 ^ 

B = ( 51 °glo( £>L ( Q ' fe '^)) ~ m obs(zi)) ^ 

C = ]T-L • (36) 

Zr « 



In the above Dl(z) = HQdi{z)/c and U{ are the errors. 
For the case of BAO data x 2 is defined as 



xi i0 (^,M) - \m^_^t (37) 
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where 



A(n° m ,a,b) = 



1 r z i dz l 2 / 3 

(38) 



Zl Jo H(z)/H 



The experimental value of A = 0.469 ± 0.017 [5] in the above is obtained from the Sloan Digital 
Sky Survey (SDSS) [5] data at z = 0.35. 
For OHD data, x 2 is defined as 



Xohd — X! 



r H th (a,b;z,i) - H ohs (zi^ " 



(39) 



In the above, <7j denote the errors and H t h(a, b; z) is obtained as (using Eqs. 8 and 31) 



H th (a,b,z) - H o {1 + 2az + {ab + a _ b)z2) ■ (40) 

The parameters a, b (and f^) are obtained by minimizing the \ 2 defined above with SNe la, OHD 
and BAO data together. As described earlier in Sect. 1 we consider five data sets of SNe la. Thus 
we have five combined data sets of (SNe la + OHD + BAO) and from the x 2 minimization five 
sets of parameters (a, b, 0,^) are obtained [14]. Needless to mention that for all the five combined 
data sets OHD and BAO data remain the same. 

4.1 Variation of u(z) with scalar fields 

The dark energy equation of state obtained from the present quintom scalar field model is denoted 
as co(z). The variation of oj(z) with scalar fields can be obtained using Eq. (30) with a chosen form 
for the scalar potentials. In the present work the form of the potentials for the fields <p and a are 
taken as 

y(0) = exp[A 1 (l + (0/M p/ ) 2 )] 

and 

V(a)=exp[\ 2 (l + (a/M pl ) 2 )} 

where M p i(= 1/8ttG) is the Planck mass and Ai and A2 are two parameters. 
From Eq. (11), we obtain for quintessence scalar field 



(m-m/h* = -f (ir^r* 21 ^ ■ (41) 



and from Eq. (18) we get for phantom field as 
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where <po and <jq are values of the fields at present epoch. With the parameters a and b obtained 
from the x 2 fit) the Hubble parameter H(z) appearing in Eqs. (41, 42) can be computed from Eq. 
(40). 

As the quantities and H a cannot directly be extracted from data, it is difficult to obtain <p(z) 
and a(z) form Eqs. (41, 42). It is also not possible to choose the boundary conditions of both the 
fields from observations. In order to circumvent this problem we propose the following formalism. 
For a varying dark energy density uj(z) with redshift z, when uo(z) > — 1, the quintessence scalar 
field (f>(z) >> cr(z) and in the case when uj(z) < —1, the phantom scalar field cr(z) >> <j)(z). In 
the former case therefore (oj(z) > —1), cr(z) can be taken to be zero for all practical purposes and 
similarly for the latter case (oj(z) < —1), 4>(z) can be taken to be zero. Thus, defining z c as the 
value of redshift z at which uj(z) = — 1, we can write, using this formalism, that in the redshift 
limit when z > z c , a(z) = and when z < z c (z > 0), 4>(z) =0. At z = z c , however, (p(z) = a(z) in 
order that the continuity of the variation of uj(z) with z is not affected. 

Under this formalism, for z > z c , the quintessence scalar field can be derived from Eq. (11) as, 

where H 2 = H\ + ^j^-p-y (neglecting H% in Eq. (29)). Also, for the case z < z c , z > 0, we get from 
Eq. (18) the phantom field as 

where H 2 = H 2 + ^^p 7 (neglecting H 2 in Eq. (29)). In the present work we are considering 
the redshift limit < z < 1.76 — the present observational reach of SNe la. In this epoch the 
universe is supposed to be dominated by dark matter and dark energy and therefore p 7 in Eq. 

(29) denotes the matter density p m . Hence the matter density can be written in terms of redshift 
z as ^rf^-pm = i?o^m(l + z ) 3 > where Ho and 0,^ are the Hubble parameter and matter density 
parameter in the present epoch respectively. 

Now as mentioned earlier in the beginning of this section, oj a ^{z) can be calculated from Eq. 

(30) with potentials exp[Ai(l + {{4> - <p c )/M p i) 2 )] and exp[A 2 (l + ((a - a )/M p i) 2 )} for different 
choices of Ai and A 2 . 

5 Results and Discussions 

Using the formalism described in Sect. 4, the x 2 tit was performed [14] for each of the five sets of 
SNe la data alongwith BAO and OHD data and the best fit values of the parameters a, b and £l° m 
obtained for all the five sets [14]. The best fit parameters for each set (alongwith BAO and OHD) 
are given in Table 1. From the obtained best fit values of the parameters the equation of state 
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MLCS+BAO+OHD 



SALT2+BAO+OHD 




ujx(z) is computed using Eq. (6) and Eq. (31). Here we mention again, that no chosen form for 
ux(z) was adopted and ux(z) was obtained in a model independent way. 



SNe la datasets 


CL 


b 


m 


(+BAO+OHD) 








HST+SNLS+ESSENCE 


1.437 


0.550 


0.268 


(+BAO+OHD) 








(Data Set I) 








SALT2(+BAO+OHD) 


1.401 


0.542 


0.272 


(Data Set II) 








MLCS(+BAO+OHD) 


1.401 


0.653 


0.296 


(Data Set III) 








UNION(+BAO+OHD) 


1.635 


0.699 


0.268 


(Data Set IV) 








UNION2(+BAO+OHD) 


1.289 


0.458 


0.272 


(Data Set V) 









Table 1. Best fit values from analysis of different data sets 

It is seen that except for Set I, for all other sets, ojx(z) extends below the limit gjx(z) = — 1 
[14]. 

As discussed earlier, the dark energy equation of state {oj(z)) is also formulated from a quintom 
scalar field model in this work. This model is based on the fact that when oj(z) > —1 its nature can 
be described by assuming a potential V((f>) for a scalar field <p responsible for the dark energy. While 
in the scenario when oj(z) < — 1 the dark energy is described by a potential V(a) for a phantom 
field a. This is described in Sections 3.3 and 4.1. Our purpose in this work is to compare dark 
energy equation of state (uj(z)) from quintom scalar field model to that (ux{z)) obtained from the 
data and check ow well the present formalism for dark energy can explain the variation of ojx(z). 

From the the best fit values of the parameters tabulated in Table 1, <1l(z) and subsequently 
Hubble parameter H(z) are computed. As discussed in Sect. 3.3, the fields <p(z) and <j{z) can 
then be calculated using Eqs. (43,44). Consequently the chosen forms of V{<p) (= exp[Ai(l + 
((4> — <p c )/Mpi) 2 )]) and V{a) (= exp[A2(l + ((<r — <7o)/M p z) 2 )]) are calculated for the particular case 
when Ai = A2 = 1. The analyses of the five sets of data [14] indicate that for Set II — Set V, we 
need to consider quintom scalar field (since cox(z) goes below —1) while from data set I only the 
quintessence field suffices. 

It is argued earlier (Sect. 4.1) that 4>(z) = <fi(z c ), when ojx(z) = —1 and (j)(z) = at other 
values of z < z c for which ujx(z) < — 1. Hence, <p{z c ) is used instead of <po in the entire calculation 
for the case z > z c . On the other hand, ao, the value of u{z) at z = 0, remains undetermined and 
hence in our calculations, we rescale a{z) as (cr(z) — ao) and show our results. 
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We use the best fit values of the parameters given in Table 1 for all the five sets of data 
and compute uj(z) formulated from quintom model following Eqs. (11,18,30). These are done for 
different values of z in the range < z < 1.76. For the present work, as stated earlier, Ai = 1 = A2 
is chosen in the expression of the potential V ((f)). They are then compared with corresponding 
ujx{z) for the five data sets, obtained from the x 2 fit an d using Eqs (6, 40) [14]. The results are 
shown in Fig. 1. 

In Fig. 1 the red coloured (solid line) plots represent the variation of ujx(z) vs z as obtained 
from the combined x 2 analyses of SNe la, BAO and OHD data [14]. The blue plots (dot-dashed 
line) in Fig. 1 are obtained from the calculation of oj(z) with the present formalism. As described 
earlier, for the first plot of Fig. 1 (corresponding to data set I), only quintessence scalar field 
describes the variation of uj(z) with the redshift z since oj(z) > —1, always, for this particular case. 
For the rest four plots the data analyses results show that ojx(z) < — 1 for some redshift values. 
Hence for comparison with oj(z), the latter is calculated for quintom scalar field involving both 
4>(z) and a(z) for those four cases in Fig. 1 (corresponding to data Set II to Set V). In these four 
plots oj(z) from our calculations are given in bi-colour plots where the blue coloured dot-dashed 
line is obtained from the quintessence field <j)(z) (oj(z) > —1) and calculations with the phantom 
field {(t(z)) are represented by only black dashed line (oj(z) < —1) of the plot. 

From Fig. 1 one readily sees that the variations of the dark energy equation of states calculated 
using the adopted forms of V(<p) and V(a) in the present quintom scalar field formalism, are 
continuous. There are no discontinuities even when uj(z) varies from <jo(z) > — 1 regime to the 
region where oj(z) < —1. More importantly, it is also evident from Fig. 1 that the dark energy 
equation of states co(z) calculated from the proposed quintom scalar field theory in the present 
work are in good agreements with cox(z) for at least three cases of data sets (namely Sets II, III 
and V) including the most recent UNION2 data set considered here. These results indicate that 
the proposed quintom scalar field formalism can be a viable model for explaining the varying dark 
energy. 
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